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In this paper the effects of time-dependent Newton constant G during inflation are studied. We 
present the formalism of curvature perturbations in an inflationary system with a time-dependent 
Newton constant. As an example we consider a toy model in which G undergoes a sudden change 
during inflation. By imposing the appropriate matching conditions the imprints of this sharp change 
in G on curvature perturbation power spectrum are studied. We show that if G increases (decreases) 
during the transition the amplitude of curvature perturbations on large scales decreases (increases) . 
In our model with a sudden change in G a continuous sinusoidal modulations on curvature power 
spectrum is induced. However, in a realistic scenario in which the change in G has some finite time 

■ scale we expect these sinusoidal modulations to be damped on short scales. The generated features 
may be used to explain the observed glitches on CMB power spectrum. This puts a bound on AG 

^ , during inflation of roughly the same order as current bounds on AG during the entire observed age 

■ of the universe. 

CO 
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i 1' I. INTRODUCTION 

o 

There are strong evidences that Newton gravitational "constant" G did not change considerably during most of the 
history of our universe. The evidences range from studies of orbits of planets to pulsar timing, evolution of stellar 
objects, the cosmic microwave background (CMB) and the effect on the abundance of light elements formed during 
the period of big bang nucleosynthesis. These studies yield a bound \G/G\ < 10~( 11-13 ) yr _1 [1, 2]. This seems to be a 
strong bound but actually it corresponds to | AG/G\ < 0.1 over the entire history of our universe, which makes it less 
dramatic. In addition, it should be noted that many other assumptions are hidden in most of these studies (like taking 
the other 'constants' to be constant). Of course, it is always possible to come up with a complicated theory of variable 
G that explains all these observations, but least deviations from the standard working model is usually preferred. 
Thus the continued lack of evidence for a variable G has been convincing for most people to accept G =constant. 
There is, however, no direct indication from any of the mentioned experiments that G has to be constant at very 
(yT) . early times, for example, during inflation. This is one of our motivations in this paper to pursue implications of the 
f^i ' change in Newton constant during the inflationary era. 

Inflation has emerged as the leading paradigm for early universe and structure formation [3J which is strongly 

■ supported by cosmological observations Q. The simplest models of inflation predict almost scale invariant, almost 
\Q I Gaussian and almost adiabatic perturbations on the CMB. In this work we would like to examine the effects of a 

■ time-dependent Newton constant in inflationary predictions. We present the formalism for curvature perturbations 
in inflationary models with a time-dependent G. As a particular example, we consider the toy model in which the 
reduced Planck mass Q = 1/8ttG undergoes a sharp change from Q = Q_ to £1 = Q + in which both fl± are constant. 

J> ' We put this change in il at an early stage of inflation so the effects of change in Planck mass is within the CMB 
observational window. As a result one expects to find local features to be imprinted on curvature perturbations from 
this sudden change in Planck mass. Indeed there has been a lot of interest in the literature to consider the effects of 
local features in inflation. These features may originate from a sudden change in slow-roll conditions, sudden change 
in the inflaton mass, particle creation during inflation, field annihilations during inflation, change in sound speed of 
perturbations or change in fluids equation of state [5143 1|. The observational motivations behind these models are to 
address the glitches in curvature power spectrum on scales associated with I ~ 20 — 40. 

The simplest implementation of a variable G is a scalar-tensor theory of gravity. In fact, there are many theoretical 
models, like string theory, supergravity and theories of extra dimensions, which involve scalar fields non-minimally 
coupled to a rank-2 tensor. Thus there is a natural place for theories of variable G within the framework of fundamental 
physics. But to maintain a constant G, and hence respect the null experiments mentioned above, authors often choose 
to change variables to another rank-2 tensor which is minimally coupled to matter. In other words, the metric that 
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describes the geometry of space-time is taken to be the one that appears in the Einstein frame. We will not do so. 
Instead we work with a model where the effective gravitational coupling, given by 51 + , has no variation during the 
standard hot big bang cosmology after inflation but experiences a change during the inflationary period. Having this 
said, if one insists on performing the field redefinition and goes to the Einstein frame with a constant G throughout, 
then the nature of the physical phenomena is the same in both frames; it is only a matter of interpretation of what 
one measures. In particular, the observers may use different clocks and rulers to measure physical processes in these 
two frames, though the final conclusion is the same for both observers. 

The effects of time-dependent Newton constant during inflation via a conformal coupling with multiple inflation 
fields to gravity was recently studied in [32| , see also [33M45j . 

The rest of this paper is organized as follows. In Section [TT] we present our setup and the background inflation 
equations. In Section IlIII we present the cosmological perturbation theory for general models with a time-dependent 
Planck mass. In Section IIVI we consider the special case in which 51 undergoes a sudden jump. By imposing the 
appropriate matching conditions we find the curvature perturbations power spectrum which will be used to put some 
rough bound on AG/G during inflation in this toy model. Discussions are given in Section W\ followed by some 
technical issues in Appendix A. 



II. INFLATION BACKGROUND 



In this section we present the inflation background. We are interested in a generic model of inflation with a 
time-dependent Planck mass with the action 

S = J d 4 xV=g(^R + C M ) , (1) 

where Cl(x a ) is a positive dimension-2 scalar function of space-time which is interpreted as the reduced Planck mass 
related to Newton constant G via 51 — l/8irG. Also R is the Ricci scalar and Cm represents the matter sector 
Lagrangian. 

The standard Einstein general theory of relativity corresponds to the case in which 51 = Mp is constant. We do not 
specify the mechanisms which make 51 dynamical. In principle, 51 may be a function of other fields in the system. As 
these fields run toward their minima a time-dependence in 51 is induced. For example, this is the situation in simple 
models of inflation with the inflaton field (f> in which 51 ~ 1 — £(f) 2 with £ being a conformal coupling (46[. 

The action contains the generic matter Lagrangian Cm- Our assumption is that the matter sector (the Standard 
Model fields) is minimally coupled to gravity with metric g a p. Of course, we can eliminate the time-dependence of 51 
by a change of metric field 



g a f3 g a f3 = -j^29a/3 , (2) 

to switch from the "Jordan frame" to the Einstein frame to ensure that G is constant in space and time. It is only 
the question of clocks and rulers one uses in these two frames and how to interpret the measurements from these two 
sets of clocks and rulers. 

We do not concentrate on a specific inflationary mechanism and the matter sector which supports inflation. All we 
need is that there is a period in which the universe has positive acceleration, d(t) > 0, where a(t) is the scale factor 
of the FRW universe with the metric 

ds 2 = g af3 dx a dx 13 = -dt 2 + a(t) 2 dx 2 . (3) 

In terms of the equation of state parameter w, we demand that during inflation w < —1/3 followed by a period of 
standard big bang cosmology with a radiation or matter dominated universe. Also we do not specify the mechanism 
in which inflation ends. 

With this general picture in mind the modified Einstein equation is 

51 G^ + (g^V 2 - V M V„)n = T^, (4) 

where T^ v — is the standard energy-momentum tensor associated with the matter action. 

A consistency condition for the above equation of motion is obtained by taking its divergence: 

VT P = -|v M 51. (5) 
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This means that Tjf is not conserved in the usual sense. Recall that the standard argument leading to V V T^ U = 
is based on applying the following mathematical fact to 5a/ = f Cm'- Given a diffcomorphism invariant functional 
SM[g,4>i] of a tensor field g and some other fields 0j, the divergence of SSmI&9ij,v vanishes if SSm/^i = for all i. 
In our case, however, Q, depends on the matter fields, so the equations of motion of those fields receive a contribution 
from J \QR as well as J Cm- However, if we define a total energy momentum tensor, T*°*, containing the mater fields 
contribution in ft, i.e. T*°* = -j=& J + Cm) /Sg a p, then T*°* is conserved in the usual sense, i.e. VT*^ = 0. 
Consider the isotropic and homogeneous FRW background with the stress energy tensor in the form 

Tp = P5 a p +(p + P)u a u , (6) 

in which p and P, respectively, are the energy density and pressure and u a is the fluid velocity four-vector normalized 
as u a u a = — 1. Assuming that Q = O(t) at the FRW background level, the background Einstein equations are 

„(*>+*) + ta _£, ( 7) 

n (zH 2 + 2H + A \ + Q + 2HQ, = -P, (8) 

2Q-+i?0 + fi =--(o + 3P). (9) 
a 3 

One can check that only two of these equations are independent. Also the consistency equation ([5]) can be written as 

p + 3H(p + P) = 3f2 (2H 2 + H+ = ^ClR. (10) 

In particular, in the standard situation in which Q = 0, Eq^. (|10[) results in the conventional energy conservation 
equation. 

In order to obtain inflation, we assume that the equation of state parameter w = P/p < —1/3 during inflation. In 
practice, we may take w ~ — 1, corresponding to a period of slow-roll inflation. We have not specified the mechanism 
in which inflation ends. It may be as in conventional scenarios of (p)reheating followed by a radiation dominated 
universe. 



III. PERTURBATIONS 



In this section we study cosmological perturbation theory in our setup with a general function n(x a ). We use the 
conventions of Ref. [47| . except that we work with 1Z C (which is the negative of their 72.) to describe the comoving 
curvature perturbations. The scalar perturbations of the metric are 

ds 2 = -(1 + 2A)dt 2 + 2ad l Bdx l dt + a 2 [(l - 2i[>)5 ij + 2d i d j E]dx i dx i , (11) 

The fluid four- vector, to linear order in perturbation theory, is = (— 1 — A, Vi) in which Vi = diV where v is the 
velocity potential. 



A. Gauge invariant perturbation equations 



The details of the perturbations of Einstein equations are given in Appendix [X] Defining the gauge invariant 
variables Newton potential $ and Bardeen potential $ 

$ = A- j t [a 2 {E - B/a)] , * = ip + Ha 2 {E - B/a) , (12) 
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the (00), (Oi), (ii) and ^ j) components of perturbed Einstein equations in the gauge invariant form, respectively, 



are 



^V 2 f + 6ffft(# + H$) + 30 (if + 2H^j + \v 2 8tt - 3H (hSQ + 5Sl) = -Sp (13) 
+ \ (a 2 Q) <I> - a( — ) = -{p + P)v (14) 



a 2 " ' \ a 

*-*-^) =0 (15) 



2£2 



>».? 

r2 



(3# 2 + 2iJ)$ + * + + + fi($ + 2* + 4H$) + 2ft$ - SQ - 2HSQ - (3H 2 + 2H)8Cl = 5P (16) 



in which the gauge invariant quantities v, Sp and 5P, corresponding respectively to perturbations in velocity potential, 
energy density and pressure, are defined via 

Tp=5 P -a 2 p[E-^j , SP = 5P-a 2 P^E-^j , v = v + a 2 (e - (17) 

Similarly, since fl(x a ) is an arbitrary function of space-time we also defined the gauge invariant perturbation in Sfl 
via 

SU = SQ - a 2 tl (e - — j . (18) 

Eq. ([15]) is very interesting. This indicates that in our setup with a time-dependent Planck mass, the Newton 
potential $ and the Bardeen potential ^ are not equal and we have gravitational anisotropy 

*-* = f. (19) 

In principle, the fact that $ ^ can have interesting observational consequences. Observationally, gravitational 
lcnsing and integrated Sachs- Wolfe effects are sensitive to $ — iff whereas galaxy peculiar velocity measurements are 
determined by the Newtonian potential $ (4^-(5l|. 

B. Separate Universe Approach 

It is also instructive to look at the separate universe app roach (53l . which is very suitable for SN formalism in 
calculating the curvature perturbation power spectrum |55l457| . In the separate universe approach the perturbations 
inside a horizon-size patch with length L ~ f/H are homogenized while a long- wavelength mode k with k > aH 
can vary smoothly over many patches (universes). In this picture the scale factor inside each patch is given by 
a(*,x) = a(i)e-^*' x ) . 

Here we generalize SN formalism in our setup to find an equation for the evolution of £ on super-horizon scales. 
For this aim, it is very useful to start with the conservation equation ([5]) contracted with the unit normal vector 
perpendicular to the surface of constant t 

= {l-A,-V i B), (20) 

which yields 

n%% = -|n^n. (21) 

After some calculation, this equation can be written as 

Sp - -(6tlR + CtSR) = -3H{Sp + SP) + 3(p + P)ip , (22) 

in which SR, the change in Ricci scalar to linear order in perturbations, is given by 

5R = -4-V 2 ($ - 2*) - &ii> - imi) - 6HA - \2(2H 2 + H)A . (23) 
or 
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To cast Eq. (|22f) into a useful gauge invariant equation for £, let us go to the uniform density gauge where Sp = 
and £ = — •0. In this gauge we have <5P = <5P na d in which <5P nac i indicates the non-adiabatic pressure perturbations. 
Considering the super- horizon modes k/aH — > 0, Eq. (|22| translates to 



£C - 3H5P nad + 3 A n 
tl 



(24) 



in which 



H 

* + (2H - — )* + + 2(i? + 2P 2 )$ 
H 



- I 0(2# + J) ) (C + *) - (A" + 2P 2 )^ . 



(25) 



In the limit where VL = Mp is constant, Eq. ([24]) reduces to the known result [53, [5J] that C = —H5P nSuC \/{p + P). 
However, in our model with a time-dependent gravitational coupling Q we have the additional contribution Aq. This 
indicates that £ does change on super-horizon scales even in the model with an adiabatic equation of state with 
<5Piad = 0. As a result, the change in f2 can be interpreted as a new source of entropy perturbations. Indeed the 
non-trivial behavior of entropy perturbations in models of multiple field inflation in which Q is a functions of theses 
fields was studied recently in |32j |. 

C. Equation for 1Z C 

In this section we provide a second order differential equation for the curvature perturbations on comoving surface 



Tl c = -ip + Hv = -* + Hv 



(26) 



which is very helpful for our analysis in upcoming sections. 

So far our equations were general. Now consider the special case of adiabatic fluid with known equation of state 
parameter w(t) = P j p. Combining Eqs. (|13[) and (| 16[) we obtain the following single second order differential equation 
for 7l c (for the details of the derivation see Appendix [5} 



e 



H 



n t 



2;„2 



T 



-jttr + ( 30 - 2f^ J Ti ck 



0. 



(27) 



Here c s is the sound speed in comoving gauge defined via SP C = c 2 s 5p c in which 5P C and 5p c are the perturbations in 
pressure and energy density in comoving gauge where v — and 
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T 



T~ 



-2fifi + 2Hrtft - AHVL 2 + 3c 2 f2 2 



= T~ 



-HQ - {5H 2 + 9H 2 c 2 s + 3H)Q - Hn(8H + 12H 2 + l2H 2 c 2 s 



ft + (3c 2 - l)Hfl + 2HVL 



= fl + 2HQ . 



and 5Q C , is defined as Sfl on comoving surfaces 



n 



5Q + riv = SQ + — (<J + Tic) ■ 
H 



(28) 

(29) 

(30) 

(31) 
(32) 

(33) 



One can check that in the limit where Q = Mp is constant, then Eq. (|2"T|) reduces to the standard result [52 



Ttck + ( -2^Wck + ^4-^ck 



H 



H 



(34) 
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IV. MODELS WITH A SUDDEN CHANGE IN NEWTON CONSTANT 



Having presented our perturbation equations in the previous section for a general time-dependent Planck mass 
0(t), here we present a toy model. Consider a model in which il goes under a sudden change at t = t c separated 
by two constant values and In this view, inflation has two stages, the first inflationary stage corresponds to 
time before phase transition, t < t c , with Q = £7_ and the second stage of inflation corresponds to time t > t c with 
il = Q + . In this picture both il± are constant and in particular f2 + = Mp in which Mp is the current Planck mass. 

In this approximation we can solve the system of equations (|13H16j) in the limits before the phase transition t < t c 
with = f2_ and after the phase transition t > t c with 51 = il + and glue the two solutions with the appropriate 
matching conditions on the background quantities and on 7Z C and 1Z C at the level of perturbations. 

First consider the background case. The background equations ([7]) and ([8]) near t — t c on each side (t = t~ and 
t = t+) are 

Ztt-H 2 _ = p- , (3Hl + 2i?_)fi_ = -p_ (35) 

and 

m+Hl = p + , (3Hl + 2H+)n+ = -p + (36) 

in which p_ = Pi^c) an d P+ = p(tt) an d so on. 

Now we consider the matching conditions on the background quantities. Our assumption is that across t = t c the 
geometry is continuous. At the background level this requires a(t) to be continuous. Integrating Eq. ([7J across t = t c 
we have J_ pdt = 3 J_ (Q.H 2 + ClH). We demand that there is no ^-function type singularity on background energy 
density, i.e., the discontinuity in p is not worse than the step function. This means that f_ pdt = 0. Furthermore, we 
expect that J_ flH 2 = since neither H nor are singular. As a result, we conclude that 

+ 

HQ. = . (37) 

This means that although (l has a 5-function type singularity at t = t c but H goes to zero at t = t c so rapidly that 
Hfl — ¥ at t — t c . This is a curious result of our assumption. This means that at t c the universe expansion comes to 
a halt and a(t) = 0. 

The second matching condition at the background level is obtained from integrating Eq. (JSJ) across the transition 
which yields 

J Pdt = 2 [i?n]± . (38) 

Here we use the notation that [X]_ = X(t+) — X(t~) for the quantity X. Eq. (|55]l indicates that P may have a 
(5-function type singularity at t = t c . 

There are different possibilities on p, P and H such that the matching conditions (|37| and (|38|) can be satisfied. One 
reasonable assumption can be that the intrinsic properties of the fluid driving inflation before and after the transition 
do not change. This means that the equation of state parameter w and the sound speed c s remain unchanged in both 
periods, i.e., 

w_ = w + = w , c s - = c s+ = c s . (39) 
Using P_ = wp- and P + = wp + in Eqs. ([35]) and (|36|) yields 

(40) 



H- _ H+ 



We can further assume that there is no jump in energy density so that p- = p+. As a result, from the continuity 
of w one also concludes that there is no jump in pressure and P_ = P+- These assumptions correspond to the case 
that there is no jump in the thermodynamical properties of the fluid driving inflation. Note that our assumption in 
taking p and P to be constant may be motivated from our staring assumption that it is the Jordan frame metric g a p 
(and not g a p as defined in Eq. ([2])) which minimally couples to the SM fields. In this view it is natural to expect the 
matter sector properties such as p and P to be continuous while the geometrical quantities such as H have non-trivial 



7 



dynamics due to the change in f2. However, at this level we do not provide any model in which these conditions are 
satisfied. In principle, one can cook up models in which these conditions may be satisfied dynamically. We will come 
back to this question briefly in Section |V] 

Having studied the matching conditions at the background level, we now consider the matching conditions for the 
perturbations. Physically, we expect both ^ and 1Z C to be continuous on the surface of transition. Geometrically, this 
corresponds to the assumption that the intrinsic and extrinsic curvatures on the three-surface at t — t c are continuous 
and 

[Kc]t=0- (41) 

Intuitively speaking, from the separate universe approach [53l. l54j|. on each patch we have a(t, x) = a{t)e-^ l ' x \ The 
continuity of the scale factor on each patch at t = t c requires ip to be continuous. Translated in a gauge invariant 
way, it is reasonable to expect VP to be continuous. Similarly one expects the curvature perturbations on comoving 
surface 1Z C = —ip + Hv to be continuous. The other matching conditions has to be imposed on 1Z C . The details of the 
matching conditions are left for Appendix [X] where it is shown that 



= . (42) 



Using the condition (j40|) this leads to 

K c+ =f3K c - , /?=J^i. (43) 




In the limit where there is no change in the gravitational coupling and ft — 1, we obtain the expected result that 
1Z C + = 1Z C - . Also note that /3 can be bigger or smaller than unity. 

Having presented the matching conditions, here we present the incoming and outgoing solutions for IZck- For this 
purpose, let us first solve the background system. Our assumption is that both periods of inflation are driven by a 
single fluid with constant equation of state parameter w and sound speed c s . To sustain inflation we require w < —1/3 
but for practical purposes we consider w ~ — 1 corresponding to a slow- roll inflation. Solving the background equations 
with (l± = for each inflationary stage one obtains 

p = p -\v) (44) 



which results in 



"('/>=<'- (l + ^-fa-'fc)) » H{rf)=U-(l + ^{r } -r 1 c)) </-//.■ U->) 



and 



-l 



a(v) = a c (l + ^(v-Vc)] , n(r,)=n + (l + ^±(r,-r, c )j r, > r, c (46) 

Here rj c represents the conformal time when ft undergoes sudden changes and a c is the value of the scale factor at 
that time a c = a{rj c ). Also % is the Hubble factor measured in conformal time H = a' /a = aH and T-L± represents 
the value of % just before and just after the phase transition: H- = HiVc) an d ^+ = H(rj£). 
The dimensionless parameter 7 is defined via 

(47) 



l + 3w 



We note that during inflation 7 < — 1 and for nearly slow- roll inflation 7 < — 1. The solution for 1Z C in each region 
of inflation is easily obtained by solving Eq. (|3~4"|) . Defining the Sasaki- Mukhanov variable u = a7Z c , Eq. (|3"4")l can be 
written in the following form 



cik z Uk = (48) 

a ' 
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in which the prime denotes the derivative with respect to the conformal time. Now using the background solution for 
a(rf) we have 



a 
a 



7(7- 1) 



7(7- 1) 



»7>17c 



(v-Vc + jHl 1 ) 



ri<Vc (V - Vc +7%- 1 ) 

Plugging this in Eq. (|48|) and imposing the Bunch-Davies vacuum for the initial conditions yields 

Tlci = C 1 x 1 {riYHP(xi(rj)) 
in which Hi (x) is the Hankel function of first kind, 

-1 



Ci 



20_a c \3k(l + w) 



(49) 



(50) 



(51) 



1 3(io -1) 

1/ = 7 — 



and 



Xi 



2(3iu+l) 



(77) = -c s fc (?7 - T} c + -IK.- 1 ) 



(52) 



(53) 



For the outgoing solution, the answer in given as a combination of the Hankel function of the first and the second 
kind 



n c2 = x 2 { v y \C 2 HP (x 2 (/?)) + D 2 H^ (x 2 (r])) 



in which C 2 and D 2 are constants of integrations and 



1 3(w-l) 
1/ = 7 — 



and 



2(3iu+l) 



x 2 (?7) = -c s fc (77 - r/ c + 7% + 1 ) 



(54) 



(55) 



(56) 



Note that in the limit where fl is continuous, we have C 2 = C\ and D 2 = 0. As we shall see, any non-trivial effect 
form change in Q is encoded in the ratio | c^ 2 1 ■ Also note that since T~L + ^ the variable x is not continuous at 
r\ c . Specifically, one can check that 



in which we have defined 



c s 7 



(57) 



(58) 



In this view, k c can be interpreted as the mode which leaves the horizon at the critical time rj c 
Now imposing the matching conditions (|4l"|) and (|43|) yields 



C2 = »7rx 2 (r? e ) 
Ci 

El 



(xM) Hl% (x 2 ( Vc )) - ( Xl (r, c )) (x 2 (t? c )) 



(59) 
(60) 



As usual we are interested in curvature perturbations, Vn, for modes which become super- horizon by the end of 
inflation when 77 — > 



(n ck K ck ,) = (2ir) 3 P n (k) S 3 (k + k') , Vn = ^Pn{k) 



k*_ 
2^' 



(61) 
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Using the asymptotic forms of the Hankel function we have 

-iT 



ft ck (?7 -> 0) 



T{v){C 2 -D 2 



(62) 



Correspondingly, the curvature perturbation power spectrum, Viz, at the end of inflation is obtained to be 

P R (i] = 0) = T? R1 (63) 
where Vtzi is the incoming curvature perturbation power spectrum and T is the transfer function 



T = 



C 2 -D 2 



In the slow-roll limit when v ~ 3/2, Vm is approximately given by 

H 2 



12(1 + w)TT 2 n 2 _c s 



(64) 



(65) 



Note that in the absence of any phase transition in which C 2 — C\ and D 2 = 0, we obtain the expected result that 
T = 1. Using the specific values of C 2 and D 2 as given in Eq. (f60|) yields 



4 P \k r 



(66) 



One can check that in the limit where ft- = fl+ and /3 = 1 we have T = 1. 

It is instructive to look at different limits of T. Consider modes which left the horizon before the time of phase 
transition, corresponding to k <C k c . Using the small argument limit of Bessel functions we obtain 



T(k < k c ) ~ 1 



/3 2 



1 



2{v-l) 



(67) 



where x = fc//c c . Interestingly, in models in which ft+ > f2_(f2_ < ft+) so /3 > l(/3 < 1), we see a shortage (excess) of 
power spectrum on large scales. This is reasonable, since for a given source of energy density, i.e. p_ = p+, we expect 
to have less gravitational interaction when the effective Planck mass increases, i.e., when + > 0_. In particular, 
when O — > oo, we expect the gravitational decoupling limit in which no metric perturbations are excited. For a view 
of T(k) see Figure [0 

Now consider the small scale limit, modes which are deep inside the horizon at the time of phase transition 
corresponding to k 3> k c . Again, using the large argument limits of Bessel functions one obtains 



T(k > k e ) ~ /3 1 



l + (/3 2 -l)cos 2 (fix- V -^-\ 



(68) 



Curiously we note that for modes deep inside the horizon, k 3> k c , T has a non-decaying sinusoidal modulations as can 
be seen in Figure [TJ Physically, one may expect that modes which are deep inside the horizon should not be affected 
by the transition in VL. The non-decaying sinusoidal modulations on small scales are a result of the assumption that 
the change in Q happens instantaneously. Realistically one expects that the change in Q takes some time scale, say 
A77. As a result for modes with frequency uj > 1/Arj, one expects T to be constant with no sin modulation. On the 
other hand for larger modes with frequency uj < 1/Arj, one expects to see localized features in T. In order to verify 
this prediction, it would be interesting to construct a model in which the change in ft has a non-zero time-scale, say 
one or few e-foldings. We would like to come back to this question in a future work. 

The spectral index of curvature perturbation power spectrum, n s , is defined via n s — 1 = dlnVu/ 'dink. In the 
model with a constant Planck mass, n, is calculated to be 



n, = 4 - 2v = 1 



6(l + u;) 
1 + 3w 



(69) 



In order to satisfy the best fit of WMAP observation n s = 0.96, we need w ~ 1.52. Now in our model with a sudden 
jump in ft the power spectrum has sinusoidal modulations and a power law ansatz for the spectral index is not a good 
one. It is an interesting exercise to check the predictions of our scenario, such as the spectral index, the amplitude of 
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tensor perturbations, etc. using the WMAP data analysis, as performed for example in (T^Ull; which predict similar 
features in power spectrum but in a different context. 

The transfer function at large k behaves differently for ft > 1 and /3 < 1. In the former case, there is a suppression 
of power: 

gi-1v < j, < pZ-2v <1 (/3 > l) (70) 

where we have made use of the fact that the observed tilt of the power spectrum is towards red, i.e., 3 — 2^ = n s — 1 < 0. 
We may ignore the non-damping oscillations of T at large k with the arguments of the preceding paragraph, but we 
expect that the first peak/trough in Figure [1] still survives. This may be employed to address the observed glitches 
in CMB power spectrum for t ~ 20 — 40. The largest error bar in this region in WMAP data is not more than one 
order of magnitude. Thus let us assume that the suppression of Viz is not smaller than a factor 1// with / > 1, e.g., 
we may take / = 2. Then it follows that for v ~ | we have j3 < y/J. In terms of Newton constant G we can write 

AC 

l-/< — <0, (71) 

where G is Newton constant measured today. This puts an 0(1) bound on the variation of G during inflation in our 
model with a sudden jump in Newton constant. Interestingly, this bound is similar to observational constraint on 
\AG/G\ < 0.1 during the observed age of the universe as discussed in the Introduction. 
In the case (3 < 1 we have an enhancement of power: 

1 < /3 3 ~ 2u < T < P 1 - 2 " (/3 < 1) . (72) 

We can now assume that the enhancement of Viz is smaller than a factor / > 1 . Then we find /? > -^j or equivalently 

0<f<l-). (73) 



V. DISCUSSIONS 



In this paper we considered the effects of change in Newton constant G or the reduced Planck mass 57 = 1/8ttG 
during inflation. The motivations for considering a time-dependent Planck mass may come from models of high energy 
physics, like string theory, in which the effective four-dimensional Planck mass is a function of extra dimension fields 
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such as the volume modulus. In the standard big bang cosmology there are tight constraints on the variations of G. 
However, there is no direct constraint on the variation of G during inflation. 

In this work we presented the formalism of cosmological perturbation theory in the models with a time-dependent 
O. In particular, we have seen that a time-varying Q causes gravitational anisotropy in which the Bardeen potential 
and the Newton potential are not equal. This can have interesting implications for the CMB and for lensing. Also, in 
the context of the separate universe approach, we have seen that a time-dependent Q effectively acts as a new source 
of entropy perturbations. As a result, £ does change on super-horizon scales even for an adiabatic inflationary fluid. 

As a specific solvable example, we considered a toy model in which Q undergoes a rapid change from f2_ to Q + 
during inflation. After imposing the appropriate matching conditions on comoving curvature perturbations 1Z C and 
its derivative 1Z C we obtained the outgoing power spectrum in terms of the incoming power spectrum. The effects 
of change in is encoded in the transfer function T(k). On physical grounds we expect that T(k) should have 
localized features centered around k — k c , where k c represents the mode which leaves the horizon at the time of phase 
transition. This is verified in our analysis. In addition, we also see a continuous sinusoidal modulation on T(k). We 
argued that this is an artifact of our simplifying assumptions that the change in Q takes place instantaneously. In 
realistic models in which the change in Q takes some finite non-zero time one expects that for very small scales T{k) 
reaches the asymptotic value unity. 

On the other hand, if one chooses to work in the Einstein frame with a constant G, then it is natural to ask what 
mechanisms create these features on power spectrum. As we discussed in Introduction and in Sec. HU in the Einstein 
frame the matter sector will have a complicated non-standard Lagrangian. In this view, in the Einstein frame one 
expects to see some sudden changes in fluid properties which causes the local features on power spectrum. In some 
ways, these sudden changes in fluid or matter sectors may be modeled via any of the mechanisms studied so far in 
literature, i.e., in 0-|3l| in the context of Einstein frame gravity. 

In our analysis so far, we have not presented a dynamical mechanism which causes the variations in Q. One 
interesting realization for creating a dynamical time-dependence in Q is to couple the gravitational system non- 
minimally to scalar fields which vary during inflation. In this view, fi can couple either to the inflaton field or to 
additional light or heavy fields present in the model. As a concrete example, consider an inflationary model with two 
fields: the inflaton field </> and the waterfall field x with the matter Lagrangian 

CM = -\{d x f-\{d4>?-v{xA). (74) 

Here the potential V is the same as in models of hybrid inflation (5|1 |59| 

1 X 1\/T 2 1 

V = 2^ + 4A" ( * 2 ~ ^T )2 + 2 9W ■ (75) 

Furthermore, suppose that 

n = n_ + ex 2 , (76) 

in which £ is a dimensionless number. 

The picture we have for the dynamics of the system is similar to waterfall dynamics in hybrid inflation. The system 
has two distinguished periods, the time before the waterfall phase transition corresponding to <f> < <j) c = M/g and 
and the time after the waterfall, <f> > <p c . The waterfall field is very heavy during inflation so it is frozen during the 
first stage of inflation and x — 0- As a result, during the first stage of inflation ft = After <f> reaches the critical 
value 4> = 4> c , the waterfall field becomes tachyonic and rapidly rolls to its global minima y = ±M/-\/A. The quantum 
fluctuations of the waterfall field S\ play crucial roles as studied in recent works a result, at the end 

of waterfall transition, on each Hubble patch \ 2 = (^X 2 ) = M 2 /X and 

n+ = 0_ + ^— . (77) 
A 

One can make the waterfall phase transition sharp enough, so the time of change in fi is reasonably short, say one 
e-folding or so. However, it is crucial to note that it is not arbitrarily sharp and in principle it takes a finite non-zero 
time scale for the waterfall phase transition to complete. This is expected to eliminate the unwanted non-decaying sin 
modulations in T(k) for small scales. Also note that, depending on model parameters [III HI]) t ne waterfall transition 
can happen either during early stages of inflation or towards the end of inflation. Having presented this dynamical 
mechanism for change in i7, it is an interesting exercise to study this model in details and see how the waterfall 
dynamics can be employed to induce a time-dependence in Planck mass. We would like to come back to this question 
in a future work. 
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Appendix A: Perturbed Equations and Matching Conditions 

Here we provide the perturbed Einstein equations in some details. The 00, Oi, ii and i ^ j components of Einstein 
equations are 



mn(ip + ha) - 30?/> - qhvla + 3H 2 sn + 3Hsn 



e r 

a 2 



sn-sie- 20^ - 2Hne 



P 



20 



sn sn 
20 ~ 20 



(Al) 
(A2) 
(A3) 



2Q{2H + 3H 2 ) + 4ffO + 20 A + (0 + 2Hfl)A + 20^' + (20 + 6HQ)ijj - (3H 2 + 2H)Sn - 2HSQ. - <S0 = SP 

(A4) 

in which Sp and SP are the perturbations in energy density and pressure, and v is the velocity scalar potential related 
to the fluid velocity four vector via u M = (—1 — A, div). Furthermore, we have defined 



9 = a 2 (E-B/a) 



(A5) 



Playing with these equations one can obtain Eqs. (|13M16I) in the gauge invariant form. 

To find the matching conditions we can find an equation relating 1Z C to as follows. To obtain this relation it is 
very helpful to go to the coming gauge where v — and TZ C = — ip. In this gauge Eqs. (|A1|) . (|A2|) and (|A4[) can be 
written as 



Sp c = 307t c - — 0<50 c - (3F0 + 4ff 2 0)<50 c + 2007^, 



„2 



<50 c + — (* + ^ c )-20* 
H 



—(2nn c - Hsn c + sn c ) - 3{h + H 2 )sn c - 3Hsn c 



-2nn c + ta = sn c - hsh c 

sp c = -6Hmz c 

in which we have defined 

T = O + 2H O , F = O + 2HVL + 2Hfl + 6i? 2 . 
Now using the definition SP C = c 2 Sp c we obtain 



(A6) 

(A7) 
(A8) 

(A9) 



- W ' 6#0-3c 2 0+^00 



+ 



T 
FH 



T 



o 







SQ C + K c - (20 - —)m 



3Hc 2 

3H + 3H 2 + ^-j^(3HQ + 4iJ 2 0) 



sn r 



H 

F 



H' 



ZHc 2 
T 







SVLc 



(A10) 



Now we look into matching conditions. We consider the model with the sudden change in 0(i) from 0_ to + 
at t = t c as presented in Sec. IIV1 In this model = before and after the transition but O has a 5-function type 
singularity at t = t c , i.e. J_ O = + — 0_. To simplify the analysis further, we also assume that SQ C — in both 

regions which is a reasonable assumption when 0± = 0. As argued in the main text, the first matching condition is 
the continuity of 1Z C across the time of transition 



[Kc£ = o . 



(All) 
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To find the matching conditions on 1Z C we use Eq. (|A10[) . Assuming that and c s are continuous as argued in the 
main text, and Cl± = 6Q C = 0, from Eq. (lAlOp we obtain 

+ 

= 0. (A12) 

Now assuming that the fluid equation of state w remains unchanged as argued in the main text, and using Eq. (|40[) 
we end up with 




(A13) 



The two matching conditions (|A11|) and (|A13j) are our two matching conditions to express the outgoing solution of 
1Z C in terms of the incoming solution. 



H ■ 
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